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Abstract Wc study role of the Hall diffusion in the 
magnetic star-disc interaction. In a simplified steady 
state configuration, the total torque is calculated in 
terms of the fastness parameter and a new term because 
of the Hall diffusion. We show the total torque reduces 
as the Hall term becomes more significant. Also, the 
critical fastness parameter (at which the total torque is 
zero) reduces because of the Hall diffusion. 
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1 Introduction 

Accretion discs are observed in many astronomical sys- 
tems from the new born stars to the compact objects 
(e.g., neutron stars) or even at the center of the galax- 
ies (e.g., Frank, King, and Raine 2002). In spite of con- 
siderable achievements in our understanding of the ac- 
cretion processes, there are many uncertainties regard- 
ing to the driving mechanisms of the turbulence in the 
disc or possible roles of the magnetic field. Although 
the standard theory of the accretion discs neglects such 
complexities (e.g., Shakura and Sunyaev 1973), it is a 
successful theory to explain parts of the observational 
features of some of the accreting systems. 

An instability related to the magnetic field and the 
rotation of the disc is known to be responsible in gener- 
ating turbulence in such systems (Balbus and Hawley 
1991). This mechanism known as magnetorotational 
instability (MRI) has been analyzed over recent years 
in detail. In addition to this role of the magnetic field, 
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the structure of the accretion discs is significantly mod- 
ified because of the magnetic effects even in the simpli- 
fied models (e.g., Lovelace, Romanova, and Newman 
1994; Shalybkov and Rudiger 2000; Shadmehri 2004; 
Shu et al. 2007). Most of the theoretical models for the 
magnetized accretion discs arc within the framework 
of the ideal MHD, in which the diversity of different 
charged particles is neglected. If ions, neutrals and elec- 
trons are considered as separate fluids, non-ideal MHD 
effects mainly due to the drift velocities between the 
charged particles will emerge (e.g., Wardle 1999; Bal- 
bus and Tcrqucm 2001). Over recent years many au- 
thors have tried to present models in order to capture 
the basic physics of non-ideal MHD effects in accretion 
discs, in particular Hall diffusion (Wardle 1999; Bal- 
bus and Terquem 2001; Rudiger and Kitchatinov 2005; 
Liverts, Mond, and Chernin 2007; Shtcmlcr, Mond, and 
Rudiger 2009). Although all these studies are showing 
significant effects of the Hall diffusion on the magne- 
tized accretion discs, the models are restricted to the 
discs around non-magnetized central object. But wc 
know that for a disc around a magnetized object like 
a neutron star, the poloidal component of the central 
magnetic field can interact with the disc and generate a 
toroidal field inside disc because of different rotational 
velocities of the central field and the disc itself (e.g., 
Ghosh and Lamb 1978). Thus, a magnetically threaded 
thin accretion disc will experience an extra torque due 
to the combined effect of the poloidal component and 
the induced toroidal component of the magnetic field. 
However, the true nature of the interaction of the stellar 
magnetic field with the disc is still under debate (Wang 
1987; Romanova ct al. 2003) and the Hall diffusion is 
neglected for simplicity. 

Following early works of magnetic star-disc interac- 
tions (e.g., Ghosh and Lamb 1978), the authors have 
been studying either the influences of the extra mag- 
netic torque on the structure of disc (e.g., Campbell 
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and Heptinstall 1998; Erkut and Alpar 2004; Matthews 
et al. 2005; Khajenabi, Shadmehri, and Dib 2008; Kluz- 
niak and Rappaport 2007) or spin evolution of the cen- 
tral star due to the exerted magnetic torque (e.g., Dai 
and Li 2006, here after DL). However, all previous stud- 
ies are restricted to the ideal MHD approach. Here, we 
study magnetic star-disc interaction, but including the 
non-ideal Hall term in the induction equation. Basic 
equations are presented in the next section. We will 
study role of the Hall term in the total torque in sec- 
tion 3. 



2 General Formulation 

The classical approach of analyzing magnetic star-disc 
interaction is done by considering thin disc approxima- 
tions for a steady state and axisymmetric accretion disc 
(Ghosh and Lamb 1978). Basic MHD equations arc 
written in the cylindrical coordinates (r, ip, z), where 
the central object with mass M is at the origin. Then, 
the equations are integrated in the vertical direction to 
obtain a simplified set of algebraic equations for the 
physical variables of the disc. We are following such 
an approach, but our analysis differs from the previous 
studies by including Hall term in the induction equation 
that will lead to a modified magnetic star-disc torque. 
Since we are interested in calculating the torque be- 
cause of the magnetic start-disc interaction, we do not 
need continuity equation or momentum equation. But 
for calculating the disc structure one should use these 
equations as well. 

First, we obtain the generated toroidal component of 
the magnetic field considering both the resistivity and 
the Hall diffusion. It will enable us to determine the 
total torque exerted on the star by the disc. Induction 
equation is 



dB 

~dt 



V x [(v x B) 
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where B = B /B and 77 is the Ohmic coefficient and Hall 
coefficient is r/n = cB/4iren Q where n Q is the number 
density of the electrons. 

In steady state, the components of induction equa- 
tion in the cylindrical coordinates (r, ip, z) are written 
as 
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Fig. 1 Geometry of the magnetic field lines in our model. 
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Fig. 2 The non-dimensional torque r /(My/ GM Ro ) versus 
fastness parameter oj (top plot). Each curve is labeled by 
the non-dimensional Hall term jn- Bottom plot shows the 
critical fastness parameter uj c (at which the total torque 
becomes zero) versus jh- We take 7 = £ = 1. 
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where B R , B v and B z are components of the magnetic 
field B. It is assumed that the magnetic diffusivity rj is 
of the same physical origin as the turbulent viscosity in 
the standard model of the accretion discs (Shakura and 
Sunyacv 1973). Thus, in analogy to the a-prescription, 
the magnetic resistivity is written as (e.g., Bisnovatyi- 
Kogan and Ruzmaikin 1976) 



n = mc s H, 



(3) 



where c s and H are the sound speed and half thickness 
of the disc, respectively. Also, r] < 1 is a dimensionless 
numerical factor. The Hall coefficient can be rewritten 
in terms of the Alfven velocity va = BJ y/4irp and the 
Hall frequency uj r = eB/cm* as 



lA. 



(4) 



where m* = p/n c is the effective ion mass and Ljj = 
«a/wh is the Hall length. 

Since we have B R < B Zl B v < B z and H << R, 
it is possible to simplify components of the induction 
equation as (e.g., Erkut and Alpar 2004) 



dB R _ dB^ v R B z | m B R d 
dz dR n n R dR 



(RB V ) = 0, 
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0. (6) 



Note that last terms of both the above equations are 
corresponding to the Hall conductivity. Equation ([5]) 



gives us 



PrB z , 
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where (3 R = —v R H/rj < 1 (Lovelace, Romanova, and 
Newman 1994), and B^ is the radial component of the 
magnetic field at the surface of the disc (see also, Shu 
et al. 2007). 

Assuming that the star is rotating with angular ve- 
locity 51, , we can integrate equation ([6]) in the vertical 
direction 
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or 
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where 7 is a parameter of order unity and ft is the ro- 
tational velocity of the disc. Here, ft K = (GM/R 3 ) 1/2 
is the Keplcrian rotational velocity. If the vertical com- 
ponent of the magnetic field threading the disc is given, 
the toroidal component is obtained using this equation. 
If we neglect the contribution of the Hall conductivity, 
we can see that equation © reduces to the classical gen- 
erated toroidal field (e.g., Ghosh and Lamb 1978). To 
evaluate B+ it is convenient to express rj^/r] in terms 
of their values at the inner edge of the disc Rq. 

In doing so, we note that the vertical component of 
the magnetic field B z due to the central neutron star 
with radius R* is 



B^B^f 
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where B* is the surface magnetic field strength. If we 
introduce nondimcnsional variable x = R/ Rq and /i = 
B*Rl as the dipole moment of the neutron star, the 
dipolar magnetic field of the central star can be written 
as 



B 

B z = — , 



(11) 



where B = fJ,/R%. Considering equations ([3]), and 
(fTTj) , the toroidal component B+ becomes 
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3 The Total Torque 

The total torque exerted on the central star because of 
the magnetic interaction with the disc is given by 



T = Tq 



' mag 1 



(14) 



where tq is the magnctospheric torque and r mag is the 
magnetic torque due to the interaction with the sur- 
rounding accretion disc. The magnetic torque r mag is 
because of the generated toroidal component B+ and 
the external vertical component B z , i.e. 



R 2 B+B z dR. 



(15) 
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But the shearing motion between the corotating magne- 
tosphere and the non-Keplerian boundary layer in the 
disc will lead to the magnetospheric torque tq. Thus, 



Ho 



-Ro-A 



i? 2 B z dR. 



(16) 



where A is the width of the boundary layer which is 
thin comparing to the inner radius Rq. To evaluate 
the magnetic torque in equation (|15|) a Keplerian rota- 
tional profile for the disc is considered. But the bound- 
ary layer does not rotate with a Keplerian profile (e.g., 
Erkut and Alpar 2004). Thus, the torque t can be 
determined by considering an approximate power law 
form for the rotational profile, 



Kq 
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Having equations (fTTj) , (|12[) and ((l"7|) . one can simply 
calculate the torque tq using equation (fl~6|) . i.e. 
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where S = A/R (5 <C 1) and u) = O*/O . For u) < 1, 
the total torque r can be written as 
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that determines relative importance of the Hall term. 
Under these conditions one can neglect the Hall term in 
the magnetized star-disc interaction if the Hall length 
is much smaller than the thickness of the disc. Other- 
wise, Hall term and its effect in reducing the magnetic 
torque should be considered as we showed in our sim- 
ple model. We also note that the Hall length given by 
equation (0| depends on the density of the electrons as 
n~ l . Thus, as the system tends to a weakly ionized 
state (i.e., density of the electrons decreases), the Hall 
length becomes larger. However, it is not easy to cal- 
culate the electron density properly because of a wide 
variety of the physical processes which are operating 
in the inner parts that result in ionization of the gas. 
For example, radiation of the central star may ionize 
the gaseous disc at the boundary layer between the disc 
and star itself. It is interesting to speculate that level 
of ionization is changed because of possible variations 
in the intensity of the radiation of the central star. If 
so, the total torque is modified because of a change in 
the number density of the electrons. Does such a varia- 
tion in the total torque modify the spin evolution of the 
central star? That is an interesting question for future 
work. But our model, at least qualitatively, shows this 
scenario is feasible. 



4 Conclusions 
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where £ = V2"fS and also we have assumed the ram 
pressure of the accreting matter at the inner radius Rq 
is equal to the magnetic pressure due to the dipolar 
field of the neutron star. So, ^ 2 /Rl = My/2GMR . If 
we neglect the Hall contribution, the above total torque 
reduces to equation (11) of DL, in which the magnetic 
star-disc interaction was studied, but without Hall dif- 



We studied magnetic interaction between a central 
magnetized star and its surrounding accretion disc con- 
sidering the Hall term. When the level of ionization is 
not high, one can expect the Hall term contributes to 
the generated toroidal field. Using a simple model, but 
illustrative, we calculated the total magnetized torque 
between the star and the disc. We found that the net 
total torque reduces when the Hall term becomes large. 
In particular, when the Hall length is a fraction of the 
disc thickness, modification to the total torque because 



'. r^n , n . -i, //TiV rPf ^he Hall term is more significant. However, we have 

b igurcl2Jshows the non-dimensional torque t/ (M y/GAl Rq)W 

r f . . . N , . not studied the structure ot the disc with this extra 

versus fastness parameter oj (top plot), Each curve is 

labeled by the non-dimensional Hall term jh and like 
DL we assume 7 = £ = 1 . The net total torque reduces 
as the Hall term becomes more dominant. On the other 
hand, we can define the critical fastness parameter cj c , 
at which the total torque is zero r = 0. Bottom plot 
of Figure [2] shows lo c versus jh- This Figure shows 
that as the Hall term becomes larger, then the critical 
fastness parameter reduces. Evidently, the importance 
of the Hall term is evaluated via nondimcnsional pa- 
rameter -fn in our model (see equation ((TBI ). If we 
assume /3r « 770 < 1 and also c s « va, then it is the 
ratio of the Hall length and the thickness of the disc 



torque which is an interesting topic for further study. 

In our study, the magnetic star-disc interaction is 
based on a direct generalization of Ghosh & Lamb 
(1978) to include the Hall term. In this approach, the 
magnetic field of the central star threads the disc at 
different radii and the field lines are closed. This model 
implies existence of an inner boundary layer extended 
to an outer transition region. However, other magnetic 
configurations have also been studied by the authors to 
analysis the magnetic star-disc interaction. For exam- 
ple, Lovelace et al. (1995) proposed a different picture 
for the magnetic star-disc interaction, in which the field 
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lines are not closed at some regions of the disc. They 
argued the magnetic field lines threading the star and 
the disc may experience a rapid inflation if the angular 
velocity of the star and that of the disc differ signifi- 
cantly. Thus, the field lines near to the star are closed, 
but far from the star the lines become open. Magneti- 
cally driven outflow or winds may launch from surface 
of the disc at the regions with open field lines and then, 
not only the mass but also the angular momentum are 
extracted from the disc by the winds. In our model, just 
the very inner part of the disc (i.e., boundary layer) and 
its magnetic interaction with the star are considered. 
In this region, the models of Lovelaec et al. (1995) and 
Ghosh & Lamb (1978) both are assuming the field lines 
are closed. However, we did not study the structure of 
the disc and only the magnetic torque analyzed. Our 
work shows that the Hall term contributes to the gener- 
ated toroidal magnetic field significantly and so, in the 
future works, one can include the Hall term and do a 
similar analysis to Lovelace et al. (1995) to study even 
the regions far from the central star and the geometry 
and structure of the winds and their back reaction on 
the disc. It is beyond the scope of the present paper. 
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